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Abstract. We consider a generalized version of the Random Energy Model in which 
the energy of each configuration is given by the sum of TV independent contributions 
("local energies") with finite variances but otherwise arbitrary statistics. Using the 
large deviation formalism, we find that the glass transition generically exists when local 
energies have a smooth distribution. In contrast, if the distribution of the local energies 
has a Dirac mass at the minimal energy (e.g., if local energies take discrete values), 
the glass transition ceases to exist if the number of energy levels grows sufficiently fast 
with system size. This shows that statistical independence of energy levels does not 
imply the existence of a glass transition. 
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1. Introduction 

In the context of condensed matter physics, the glass transition is a rather generic 
phenomenon through which the state of a system becomes partially frozen below a 
threshold temperature [H [2j [3] . One of the simplest models exhibiting a glass transition 
is the Random Energy Model (REM) introduced by Derrida [I] , in which the energies of 
microscopic configurations are independent and identically distributed random variables 
drawn from a distribution V(E), chosen to be Gaussian in the original version of 
the model. The REM provides a simple illustration of the so-called 'one-step replica 
symmetry breaking' scenario of the glass transition [5], which is known to hold in more 
sophisticated mean- field models P, [7]. Several generalizations of the REM have been 
proposed, mostly to incorporate correlations in simple ways [H [TUHH]. In addition, 
interesting connections of the REM with probabilistic issues such as the convergence 
properties of sums [12] and extreme values [131 QUI EH] of random variables, as well as 
with signal processing issues such as moment estimation [15], have been pointed out. 
Potential connections with string theories have even been recently outlined [T6j . 
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In this note, we explore the question whether uncorrected random energy levels 
are enough to generate a glass transition. To this aim, we consider a generalized version 
of the REM offering some freedom in the energy distribution as well as in the scaling 
of the number of configuration with system size. The paper is organized as follows. 
The model is introduced in Sect. El and the necessary framework to analyze the glass 
transition is presented in Sect. |3j Then in Sect. HI which constitutes the core of this 
paper, we derive necessary conditions for the absence of glass transition, and show 
that these conditions are also sufficient. In addition, we determine the behavior of the 
glass transition temperature close to onset. Finally, Sect. [5] summarizes the results, and 
discusses some analogies with related problems. 

2. Model 

2.1. Definition 

We consider a disordered system having a number M of microscopic configurations, 
labeled by index k — 1, . . . , M; to each configuration is associated a quenched random 
energy E^. The energies are assumed to be independent and identically distributed. 
At variance with the standard REM, we do not directly specify the distribution V(E) 
from which the energies Ej, are drawn, but we rather assume that E\. is given by a sum 
of N individual contributions, 

N 

E k = J2vk,i, k = l,...,M, (1) 

i=l 

where the terms 77^ are independent and identically distributed random variables with 
a finite variance distribution p{rj). A standard assumption is that the number M of 
configurations scales as 

M ~ e aN (a > 0). (2) 

For the standard REM, a = In 2. Note that this model may be considered as a simplified 
version of the directed polymer problem in a random media jTTj, neglecting correlations 
between local energies on different paths. 

One can interpret Eq. flTJ as the decomposition of the total energy over the N 
degrees of freedom of the system, with the strong assumption that local contributions 
associated to different microscopic configurations k are statistically independent. Such 
an approach provides an alternative interpretation (besides the standard one in terms 
of p-spin model in the limit p — > 00 [4]) of the Gaussian energy distribution of 
variance proportional to N used in the standard REM jl] (though more general energy 
distributions have also been considered [5]). But it also gives many ways to depart 
from the Gaussian distribution, in the sense that equilibrium at finite temperature is 
dominated by energy values in the lower tail of the energy distribution, far from the 
maximum of the distribution around which the Gaussian approximation holds. As we 
shall see below, the key ingredient of our generalized version of REM is the introduction 
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of both the arbitrary distribution p(rj) and the free parameter a, which leads to a richer 
behavior than in the standard REM. 



2.2. Large deviation function 

To characterize the distribution V[E) of the energies Ey. from the distribution p(rj), we 
use large deviation theory [TBI [19]. Let us define the energy density e = E/N, and its 
distribution P(e) = NV(Ne). The Gartner- Ellis theorem [181 EH] implies that, if the 
cumulant generating function 

X(q) = In (e™> (3) 

is defined and differentiable on the real axis, then the distribution P(e) is given by 

P(e) « e"^ (£) (4) 

where is the Legendre transform of A, 0(e) = max g {eg— A(g)}. Note that the properties 
of the Legendre transform imply that 0(e) is a convex function. Similarly, since for all 
e, (eg — X(q))\ q= o = 0, we have 0(e) > 0. Moreover, the lower bound is attained at 
e = A'(0): 0(A'(O)) = 0. In the present statistical physics context, it is however more 
natural to use the function 

M/3)=hi<e-^)=AR3) (5) 

instead of A(g), and we shall thus use fi(/3) in the following, together with the relation 

0(e) = -min{e/3 + /i(/3)}. (6) 



3. Analysis of the glass transition 

In the present section, we briefly set up the general framework allowing the glass 
transition to be studied. This framework mainly relies on the existence of a finite 
size cutoff in the density of state jl]. 



3.1. Finite size cutoff 

Considering a given sample with a finite number of configurations, the energies Ej~ 
are necessarily confined to a finite subdomain of the support of V(E). In the low 
temperature regime, the lower bound of this domain is known to play an important role 
[I] . Although the boundaries of this domain are also random variables, it is nevertheless 
possible to define a sharp lower boundary on the energy density in the large M limit. 
Defining the cumulative F(e) = f_ de' -P(e'), we consider the probability 

P(e 1 ,...,e A/ >e) = [l-F(e)] M (7) 

that all energy densities e^, k — 1, . . . , M, are larger than a given value e. We denote as 
eo the value of e for which P(e) is maximum (i.e., 0(e o ) = 0). Note that 0(e), being a 
convex function, necessarily decreases for e < eo- Using a standard saddle-point method 
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and neglecting non-exponential prefactors, we have for e < eo, F(e) ~ e N ^ so that 
[1 - F(e)] M « e" MF(€) . Eq. (JT)) can then be rewritten as 

lnP( ei ,...,e M > e) « -MF(e) = _ e ^-^)), (8) 

using InM = a A" [see Eq. (j2J)]. Let us define < e such that 

0(et) = a. (9) 

Eq. (IH]) shows that P (ei, . . . , tu > e) exhibits a sharp crossover at e = e* . For e < e\ 
a — 0(e) < 0, and P(ei,...,ej^ > e) — > 1 when iV — )■ oo, so that there is with 
probability one no energy levels below e. In contrast, for e > e* , a — 0(e) > 0, so 
that P(ei, . . . , 6m > e) — ► when AT — >■ oo, which means that with probability one, 
there are energy levels lower than e. In the limit N — > oo, the value e^ thus corresponds 
essentially to a border below which there are no more energy levels. In other words, the 
value e^ can be considered as the ground state energy density in almost all samples 0. 

3.2. Free energy and glass transition temperature 

To evaluate the disorder averaged free energy (F) = —(3~ 1 (lnZ), where Z is the 
partition function Z = J2k=i e ~ l3Ek an< ^ P = V^ 1 ^ ne inverse temperature, a usual 
method (beyond the replica trick [5]) is to determine the typical value Z typ (rather than 
the averaged one) of the partition function, yielding (F) w — TlnZ typ . Z typ is evaluated 
taking into account the threshold e\ and approximating the density of states n(e) by 
the disorder averaged one (n(e)) = M P(e) in the energy range e > e\ where (n(e)) is 
large. One finds, using the large deviation form Eq. (j4j), 

e Ns ^de, g(e) = a-<f>(e)-Pe. (10) 

For large N, the partition function Eq. (TTOT) can be evaluated through a saddle-point 
approximation. Equating the derivative g'(e) to zero yields 

0'(e m ) = -(3. (11) 

As g"(e) < 0, g(e) decreases for e > e m . If e m > e\ the saddle-point evaluation yields 
Ztyp ~ e Ar ^ a_ ^ ( - em - )_/3em - ) . In the opposite case e m < e\ the global maximum of g(e) is no 
longer relevant as it falls outside the integration interval. The maximum of g(e) over 
the interval [e^, oo) is then g(^), leading to Z typ w e N(a-<j>{e^)-p^) _ ^Qj-^gj- between 
these two regimes, e m = e', defines the glass transition temperature T g = P^ 1 through 
the implicit relation 

e m (/3 g ) = et. (12) 

As 0(e) is a convex function, 0'(e) is an increasing function of e. Hence from Eq. (ITT]) . 
e m is a decreasing function of j3. The glassy regime e m < e^ thus corresponds to /3 > (3 g , 

| Note that similarly, an upper bound also exists for the energy levels, but we do not take it into 
account as it plays no role in the thermodynamic properties. 
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or equivalently to T < T g . Altogether, the free energy per degree of freedom / = (F) /N 
reads for N — > oo, using Eq. Q, 

/(/3) = ^0(e m (/3)) + e m (/3)-^, < S , (13) 

/OS) = e*, > (14) 
Inverting the Legendre transform Eq. (151) , one finds 

H{0) = -/3e m (/?) - 0(e m (/3)) (15) 

so that the free energy can be rewritten for /3 < /3 g as /(/?) = — i /i(/3) — § . Introducing 
the function 

C(/3) = /V(/3)-M/3), (16) 
the entropy s = —df/dT = (3 2 df/d(3 then reads 

««_{.-«». j<j (17) 

Note that C'(/^0 = Pv>"{P) and that is a convex function with yu(0) = 0, so that 
(((3) is an increasing function of f3, starting from £(0) = 0. Accordingly, the entropy s 
is a decreasing function of (3 for /3 < /3 g . We thus recover in this general framework the 
standard interpretation of the glass transition in terms of a vanishing entropy per degree 
of freedom, meaning that in the low temperature phase, the probability distribution 
concentrates on a few microscopic configurations [I]. Eq. (JT7J) provides us with an 
alternative characterization of the transition temperature /3 g 

C()9g) = «■ (18) 
From Eq. (fl8l) . one sees that the glass transition exists if the function C(P)> which is 
defined for all /3 > 0, reaches the value a for some finite inverse temperature In the 
standard REM, which is recovered by choosing for p(rj) a centered Gaussian distribution 
of variance J 2 /2, one has ((f3) = J 2 (3 2 /4, so that (((3) can reach any value a, implying 
the existence of the glass transition. More precisely, as shown in the left panel of Fig. [fl 
we have 

= \/^- (19) 

However, in the present more general setting of an arbitrary p(rj), ({(3) may be bounded 
and the glass transition may not exist, as we shall see in the next section. 



4. Conditions for the existence or absence of glass transition 

In this section, we study the asymptotic behavior of (((3) for /3 — > oo, to see if £(/3) 
converges to a finite limit, or diverges. If C{P) diverges when (3 — > oo, it will necessarily 
cross (assuming continuity) the value a for some finite (3 g . In constrast, if C{P) converges 
to a finite limit C^, Eq. (fl8l) has a solution only if a < (c^. Hence a glass transition 
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Figure 1. Sketch of the phase diagram of the REM in the (a, /3) plane. Left: case when 
the function £(/3) diverges when (5 goes to infinity. The line corresponds to the glass 
transition (inverse) temperature /3 g (a), illustrated here on the standard REM with a 
Gaussian energy distribution -see Eq. (|T9l) . Right: same diagram in the case when the 
function C(/3) has a finite limit a c . The vertical dashed line a = a c corresponds to the 
asymptote of /3 g (a). For a > a c , no glass transition occurs. 



exists only if a < a c = Coo as illustrated on the right panel of Fig. [TJ The question is 
then to know for which form of the distribution p(rj) -see Eq. (jSJ)- the function can 
have a finite limit a c . 

4-1. Necessary conditions for the absence of glass transition 

In this section, we wish to derive some necessary conditions for the absence of glass 
transition. We thus start by assuming that C{P) has a finite limit a c and we explore the 
implications of this assumption on the distribution p(rj) of the local energies. 

4-1.1. Asymptotic behavior of [//(ft) and p((3) 

In order to derive the behavior of p(/3) from the hypothesis that C{P) admits a finite 
limit, it is useful to express /i((3) as a functional of £(/?) by solving Eq. (TIB]) as a 
differential equation in /i, for a given function (. Using classical ordinary differential 
equation method, the following result is obtained 



where p!^ is a constant. One should note that the integral in Eq. (f2D|) is always well- 
defined if £ admits a finite limit. 

The next step is to use this integral form of p to show that p!^ is rightfully the 
limit of p'(P) when (3 — > +oo. Differentiating Eq. ( 12 Op yields 




(20) 




(21) 




(22) 
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Note that from the convexity of p((3), p'{f3) is an increasing function, so that 

/A/3)-/4<0, (23) 
a property which will prove useful later on. We shall show that p((3) has a linear 



asymptote, which is actually not obvious from Eq. ( 122]) Eq. ( 120]) can be rewritten as 
KP) - /Voo = ~P J ^rdt. (24) 
Since (((3) admits a finite upper bound K, we have 

/.+oo ^ 

A*0») - /Voo > -P J fdt = -K. (25) 

From Eq. fl23|) . yu(/3) — is a decreasing function. Therefore, the existence of the 
lower bound derived in Eq. ( 125]) implies that p((3) — fip'^ has a finite limit, namely 

hm / i(/3)-/3/i' 00 = /i 00 . (26) 

/3->+oo 

In other words, /i(/3) has a linear asymptote with slope p'^. Furthermore, one necessarily 
has fioo < since p((3) — Pp!^ is a decreasing function starting from the value at = 0. 

4-1-2. Support of p(rj) 

We shall now focus on the case p'^ = 0, as the case of an arbitrary value p'^ can be 
obtained from the case p!^ = by a shift of the variable 77. We shall first specify the 
support of p{rf). One of the properties of the cumulant generating function A introduced 
in Eq. ([3]) is that the image of 1R by A' is the support of the probability density function 

{A'(/3), f3 e M} = { V , p( V ) > 0, 7] e R}. (27) 

Given that p(/3) = A(— 0) [see Eq. (|3j)], the lower bound of the support of p{rf) is directly 
related to p'^ 

7] min = M{7), p(7j) > 0} = Voo- (28) 

Notably, p'^ = implies that only the positive values of rj have a non-zero probability. 
Taking into account the fact that the support of p(7j) is [0, 00), p((3) can be rewritten 

as 

POO 

p(/3) = In / p(r])e-^d7]. (29) 



n 



§ For instance, a function of the type f(x) = ax + b In x has a finite derivative equal to a when x 
but has no linear asymptote. 
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4-1.3. Characterization of p{rj) in the neighborhood of r) m i n 

Eq. (129]) implies that the asymptotic behavior of /i(/3) (for — > oo) is directly linked 
to the behavior of p(r/) in the neighborhood of ?7 m i n = — fi'^. Moreover, as stated 
by Eq. (1261) . for jj!^ = 0, /i(/3) admits a finite limit /ioo- We shall now explore the 
consequences on p(rj) of the existence of this finite limit /ioo- 

If /ioo = 0, the property /i(0) = and the monotonicity of imply that /x(/3) = 0. 
The only possibility for p(rj) is the degenerate distribution p{rj) = 5(rj). This is obviously 
not a situation of interest, since all configurations of the model would have the same 
energy, equal to zero. 

If /ioo < 0, p(rj) has to satisfy 

lim / p(n)e-^dr] = e Moc , < e Mo ° < 1. (30) 



The integral in Eq. ( |30j) can be rewritten as, assuming /3 has been made dimensionless, 

> r^/VP roc 

p{r))e-^dn= \ p^e'^drj + / p(n)e~^dn. (31) 
We first note that the second integral in the r.h.s. of Eq. ( 13TT) can be bounded as 

/*oo /*oo 

/ p{n)e- pri dn < / p(r/)rfr/ < e"^, (32) 

and thus goes to zero when /3 — > oo for any distribution p(r/). We now focus on the 
behavior of the first integral in the r.h.s. of Eq. ( 13~1|) . If p(n) does not contain a Dirac 
mass (i.e., a Dirac delta) at n = 0, this integral goes to zero when /3 —> oo, while in the 
presence of a Dirac mass at n = 0, the integral takes a finite limit when (3 — > oo. As a 
result, Eq. f )30|) . which is a consequence of the assumption that C(/3) has a finite limit 
when /3 — > oo, can only be satisfied if p(n) contains a Dirac mass at n = (we recall 
that n = is the lowest accessible value of 77). In this latter natural form for 

p{rf) then consists of the following mixture: 

p( V ) = D5( V ) + (l-D)c( V ) (33) 

where < D < 1 and c(n) is a probability density function with support [0, 00) and 
no Dirac mass at n = 0. For instance, c(n) can be a regular distribution (which may 
have an integrable divergence at n = 0) or a sum of Dirac masses, meaning that n takes 
discrete values. 

Finally, as mentioned above, the case fi'^ 7^ can be easily obtained through a shift 
of the variable 77, so that the same results hold in full generality. The generalization of 
the distribution given in Eq. ( 1331) reads 

p(r ] ) = D5(r ] + f i' O0 ) + (1 - D) c( V ) (34) 

with the support of c(n) limited to [—(jf^, 00). 

In the following subsection, we check that the distribution Eq. ( |33l . obtained 
through necessary conditions, indeed leads to a finite a c , and thus to the possibility 
of the absence of the glass transition. 
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4-2. Distributions with a discrete mass at the minimal energy 

We now wish to show that if one starts from the distribution p(rj) given in Eq. ( )33|) . the 
resulting function (((3) converges to a finite limit a c . This result is not obvious from 
the previous subsection, where we used necessary conditions only, and did not study the 
behavior of the term 0fi'(0) in ((0)- Starting from Eq. ( 1331) . can be expressed as 

fi(0)=\n(D + (l-D)I(0)), (35) 

where we have defined 



/(/?) = / c (rj) e-^drj. (36) 



The derivative n'(0) is then given by 

»>(0) ~ ^- D)I '^ (37) 

M/ * J ~ D + (i-D)i(py [ n 

As shown in Eq. ( 13T|) . an integral of the form of I ((3) converges to zero when — > oo 
for any distribution c(r]) which has no Dirac delta at rj — 0. We thus have fi'(0) ~ 
(1 — D)I'(0)/D for —> oo. The evaluation of 0fi'(0), which is needed to compute 
COS), thus boils down to that of 01' (/3). 

The behavior of 01' (0) can be evaluated by performing an integration by part, 
yielding 

poo 1 

0I'(0) = - ^(vc^e^drj. (38) 
Jo dv 

Following the same arguments as for 1(0), one can then show that 01' (0) goes to zero 
when — >■ oo, and so does 0n'(0). Hence from Eq. (fl6|) . ((0) converges to the finite 
limit a c = — In D > when — > oo. 

In conclusion, for a distribution p(rj) of the form Eq. fl33|) . we have shown that the 
glass transition disappears for a > a c , with 

a c = -\nD, (39) 

as sketched in Fig. [2j 

4-3. Behavior of the glass transition close to the onset threshold 

We explore some interesting consequences of the above results. For a < a c , the glass 
transition exists, but the temperature range of the glassy phase is expected to shrink 
when a — > ct c . Let us make the argument quantitative. The glass transition temperature 
is determined from the relation (,{0g) = a [ see Eq. (1181) ]. For a close to (and smaller 
than) a c , g is thus determined by the asymptotic behavior of ((0) for large 0. 

To make concrete calculations, we assume that 0(77) ~ cqti u ~ 1 when 77 — >- [y > 0). 
Using the change of variables u = 0rj in the integral defining 1(0), see Eq. (136]) . one 
finds for — > 00 that 

^4f*G0"~ WE ^' (40) 
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Figure 2. Phase diagram in the (a,D) plane, for the distribution p{rj) defined in 
Eq. The line a c = — \nD separates regions where the glass transition exists 

(a < a c ) from regions where no glass transition occurs. 



can be computed in the same way, simply replacing 0(77) by —770(77), which amounts 
to replacing v by u+ 1 and Cq by — Cq. One then finds ~ —T(u + 1)cq/ f3 u+1 . Taking 
into account Eqs. ( )35l) . ( 1371) and ( )40l) . we obtain 

-ln -CW~ (1 -° )( ;; 1)rWC ° . (41) 

Using C(Pg) = «j we get 

-InO-q-t 1 - )'^ 1 ^. (42) 

At this stage, two different viewpoints can be adopted, namely either considering a 
as the control parameter for a fixed D, or considering D as the control parameter for 
a fixed a. We first fix D, and use a as control parameter. In this case, the glass 
transition occurs for a = a c = — InD. From Eq. (142|) . we obtain for the glass transition 
temperature T g = (3' 1 

T g ~ A (a c - a) 1/v , a a~ , (43) 

with A = [(1 - D){v + l)T(u)c /D]- l / u . 

Alternatively, using D as control parameter for a fixed a, the transition exists only 
for D < D c = e~ a . Eq. flU then leads to 

T g ^ A (D c — D) 1 /", D^D~, (44) 

with A=[(l- e~ a ){y + l)r(^)c ]- 1 / iy . 

5. Discussion 



In summary, we have shown that by tayloring the distribution 79(77) of local energies 
(and thus modifying the global energy distribution V(E)), the glass transition can be 
avoided in some parameter regimes, so that the presence of uncorrelated random energy 
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levels is not a sufficient condition for the emergence of a glass transition. Reversing 
the perspective, one could also interpret the onset of a glass transition when varying 
either a or the weight D of the Dirac mass, as a kind of critical phenomenon, the order 
parameter of this transition being the glass transition itself. We have seen in Sect. 14.31 
that the glass transition temperature indeed behaves as a power law close to threshold. 
However, the corresponding exponent is non-universal, as it depends on the behavior of 
the distribution p(rj) close to the lower bound of its support. 

These results further suggest an interesting analogy with another type of critical 
phenomenon. We have seen that if we build the random energies of the model by 
adding up random positive terms r)k,i, such that each term is either equal to zero with 
probability D, or drawn from a continuous distribution c(rj) with probability 1 — D, the 
glass transition disappears if the fraction of zero terms exceeds some threshold D c . This 
situation is reminiscent of the dilute Ising model, where the coupling constants between 
neighboring sites are randomly set to zero with a given probability. Above a critical 
fraction of zero couplings, the transition disappears [201 EH E2] . The mechanism at play 
in the dilute Ising model is however different, as it is related to the percolation of the 
bonds with nonzero couplings, and thus has a geometric interpretation. 

Besides, it is interesting to study the properties of the large deviation function 
0(e), which is related to the microcanonical entropy s m (e) through s m (e) = a — (f>(e) 
[see Eq. (110j) ]. Using Eq. (115j) and the properties of the Legendre transform, one can 
show that 0(e m (/3)) = C(/^)- Denoting as 77 min the lower bound of the support of the 
distribution p{rf), one has e m (/3) — > r] m i n when (3 — > oo (note that ?7 m i n may be equal to 
— oo). Hence if (((3) — >• oo for (3 — > oo, 0(e) also diverges for e — > r] min . In this case, the 
equation 0(e^) = a always has a solution > r/ m in, assuming that 0(e) is continuous. 
Furthermore, as 0(e) is regular at e = e\ its derivative is finite, and so does the glass 
transition temperature /3 g = 0(e^) -see Eqs. ffTTI) and (fT2l) . 

In constrast, if ((/3) converges to a finite limit, 0(e) also has a finite limit, equal 
to a c , for e — > ?7 m i n . In this case, as 0'(e m ) = —/3 (see Eq. (fTTj) ). the large deviation 
function necessarily has an infinite negative slope for e — > rj m i n . Hence the absence of 
glass transition in this case does not have the same origin as in non- disordered systems. 
For such systems, there is no cut-off e\ and energies down to rj m \ n can be explored. In 
the presence of disorder, the glass transition can be avoided only by tayloring 0(e) so 
that the equation determining the cut-off has no solution > ?7 m i n . 

Let us finally mention that the glass transition in the REM also has applications in 
other fields, like the empirical estimation of moments in statistical signal processing. It 
has been recently emphasized that the so-called "linearization effect" in multifractal 
analysis, occuring when empirically determined moments Sm(q) = ^~ 1 X]f=i a; ? 
significantly depart from the theoretical ones (x q ), can be interpreted as an analog of the 
glass transition in the REM [T5l 1231 1241 125] . the inverse temperature (3 being mapped 
onto the moment order q. Similar effects also occur even when considering uncorrelated 
signals, as long as the marginal distribution of the signal is sufficiently broad, but with 
finite moments (for instance a lognormal distribution) [26]. Along this line of thought, 
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the present version of the REM can be seen as the analog of the moment estimator of 
variables built as products of a large number of independent random variables (which 
is a possible way to build broadly distributed variables with finite moments). The 
present study then shows that in most cases, this moment estimator will also present a 
linearization effect, and depart from the theoretical moments for q above some threshold 
q*. However, if the underlying distribution of the variables presents a Dirac mass at 
its upper bound, this Dirac mass generates a linear branch in the logarithm of the 
theoretical moments which can mask the empirical linearization effect. 
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